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Abstract 



We consider a Josepiison junction system installed with a finite length inhomogeneity, either of 
microresistor or of microresonator type. The system can be modelled by a sine- Gordon equation 
with a piecewise-constant function to represent the varying Josephson tunneling critical current. 
The existence of pinned fluxons depends on the length of the inhomogeneity, the variation in the 
Josephson tunneling critical current and the applied bias current. We establish that a system 
may either not be able to sustain a pinned fiuxon, or - for instance by varying the length of the 
inhomogeneity - may exhibit various different types of pinned fluxons. Our stability analysis shows 
that changes of stability can only occur at critical points of the length of the inhomogeneity as a 
function of the (Hamiltonian) energy density inside the inhomogeneity - a relation we determine 
explicitly. In combination with continuation arguments and Sturm-Liouville theory, we determine 
the stability of all constructed pinned fluxons. It follows that if a given system is able to sustain at 
least one pinned fiuxon, there is exactly one stable pinned fiuxon, i.e. the system selects one unique 
stable pinned configuration. Moreover, it is shown that both for microresistors and microresonators 
this stable pinned configuration may be non-monotonic - something which is not possible in the 
homogeneous case. Finally, it is shown that results in the literature on localised inhomogeneities 
can be recovered as limits of our results on microresonators. 

Keywords: Josephson junction, inhomogeneous sine-Gordon equation, pinned fiuxon, stability . 
AMS subject classifications: 34D35, 35Q53, 37K50 . 

1 Introduction 

In this paper we consider a sine-Gordon-type equation describing the gauge invariant phase difference 
of a long Josephson junction 



where x and t are the spatial and temporal variable respectively; (/){x, t) is the Josephson phase difference 
of the junction; a > is the damping coefficient due to normal electron flow across the junction; and 
7 is the applied bias current. The parameter D represents the Josephson tunneling critical current, 
which can vary as a function of the spatial variable. 
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When D is constant (without loss of generahty, we can take D — 1) and there is no imposed current 
and dissipation, i.e., 7 = a = 0, the system (1) is completely integrable [1] and has a family of travelling 
kink solutions of the form 

(j){x^t) = 00 I — 1^=^^ I , with (/>o(0 = 4arctan(e^) for any \v\ < 1. (2) 
V V 1 — '^^ / 

In the study of Josephson junctions, this kink represents a fluxon, i.e. a magnetic field with one 
flux quantum $0 - 2.07 x 10"^^ Wb. If there is a small induced current and dissipation but no 
inhomogeneity, then there is a unique travelling fluxon whose wave speed in lowest order is given by 
V = . ^ — and no stationary fluxons exist, see, e.g., [111. 

^16(a/7)2+7r2 

It was first suggested and shown in [22] that if the critical current D is locally perturbed, stationary 
fluxons can exist even if an imposed current is present (7 7^ 0) and that a traveling fluxon (2) can be 
pinned by the inhomogeneity. This phenomenon is of interests from physical and fundamental point 
of view because such an inhomogeneity could be present in experiments due to the nonuniformity in 
the width of the transmission Josephson junction line (see, e.g., [2, 28]) or in the thickness of the oxide 
barrier between the superconductors forming the junction (see, e.g., [29, 32]). About a decade after 
the first analysis of this phenomenon, it is shown in [17] that the interaction between a soliton and 
an inhomogeneity can be non-trivial, i.e. an attractive impurity, which is supposed to pin an incoming 
fluxon, could totally reflect the soliton provided that there is no damping in the system. Recently it is 
proven that the final state at which a soliton exits a collision depends in a complicated fractal way on 
the incoming velocity [12]. 

So far almost all of the analytical and theoretical work describes the local inhomogeneity by a 
delta-like function [12, 16, 17, 22]. Yet, the length of an inhomogeneity in real experiments is varying 
from 0.5Aj [29] to 5Aj [2, 28], with Aj being the Josephson penetration depth. Therefore, such 
inhomogeneities are not well described by delta- functions. Kivshar et al. [15] have considered the 
time- dependent dynamics of a Josephson fluxon in the presence of this more realistic setup, i.e. fluxon 
scatterings that take into account the finite size of the defect, within the framework of a perturbation 
theory, i.e., when a, 7 are small and D ^ 1. Piette and Zakrzewski [26] recently studied the scattering 
of the fluxon on a finite inhomogeneity, extending [17, 12] for finite length defects in the case when 
neither applied bias current nor dissipation is present. 

In this paper, we consider the problem of a long Josephson junction with a finite-length inhomogene- 
ity and give a systematic analysis of the existence and stability of stationary fluxons, i.e., the pinned 
fluxons, as they lie at the heart of the interaction of the travelling fluxons with the inhomogeneity. The 
Josephson tunneling critical current (denoted by in (1)) is a function of space and is modelled by 
the step-function 

D{x;L,d)^h; (3) 
I 1, \x\ > L. 

with d > 0. The inhomogeneity, as modelled by (3), can be fabricated experimentally with a high 
controllability and precision, such that the strength and the length of the defect d and 2L can be 
made as one wishes (see [33, 34] and references therein for reviews of the experimental setups). When 
the parameter d is greater or less than one, the inhomogeneity is called a microresonator respectively 
microresistor. They can be thought of as a locally thinned respectively thickened junction, which 
provide less respectively more resistance for the Josephson supercurrent to go across the junction 
barrier. Note that as (1) without inhomogeneity is translationally invariant, it does not matter where 
the inhomogeneity is placed. The existence and stability problem of pinned fluxons in finite Josephson 
junctions with inhomogeneity (3) has been considered numerically by Boyadjiev et al. [3, 6, 7]. Here 
we consider an infinitely long Josephson junction with inhomogeneity (3) and provide a full analytical 
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study of the existence and stability of pinned fluxons, using dynamical systems techniques, Hamiltonian 
systems ideas, and Sturm-Liouville theory. 

For the existence of the pinned fluxons we observe that, as I? = 1 for |x| large, it follows immediately 
that the asymptotic fixed points of (1) are given by siiKj) = 7, and the temporally stable stationary 
uniform solutions are = arcsin7 + 2k7r. Bu definition, a pinned fluxon is a stationary solution of (1), 
which connects arcsin7 and arcsin7 + 27r. Hence a pinned fluxon is a solution of the boundary value 
problem 

(l)xx - D{x; L, d) sin + 7 = 0; 

(4) 

lim (f){x) — arcsin7 + 27r and lim (j){x) — arcsin7. 

First we observe that pinned fluxons can only exist for bounded values of the applied bias current, 
I7I < 1 (where this upperbound is directly related to our choice to set D = 1 outside the defect). 
Moreover, there are symmetries in this system. If (f){x) is a pinned fluxon connecting arcsin7 (at 
X —00) and arcsin7 + 27r (at x +oc), then (j)(—x) is a solution as well, connecting arcsin7 + 27r 
{x —00) and arcsin7 (x +oc). So the second solution is a pinned anti- fluxon. The symmetry 
implies that we can focus on pinned fluxons and all results for pinned anti-fluxons follow by using the 
symmetry x —x. Another important symmetry is 

0(x) 27V — (j){—x) and 7 —7. 

Thus if (j){x) is a pinned fluxon with bias current 7, then 2ti — (j){—x) is a pinned fluxon with bias 
current —7. This means that we can restrict to a bias current < 7 < 1 and the case —1 < 7 < 
follows from the symmetry above. 

Furthermore, the differential equation in (4) is a (non-autonomous) Hamiltonian ODE with Hamil- 
tonian 

i7 = -p^ — i?(x; L, (i)(l — cos(/)) + 70, where p — (l)x. (5) 

The non-autonomous term has the form of a step function, which implies that on each individual 
interval (— oc, — i^), (— i), and (L, oc) the Hamiltonian is fixed, though the value of the Hamiltonian 
will vary from interval to interval. Therefore the solutions of (5) can be found via a phase plane 
analysis, consisting of combinations of the phase portraits for the system with D — \ and D — see 
also [30] for a similar approach to get existence of vr-kinks. In the phase plane analysis, the length 
of the inhomogeneity (2L) is treated as a parameter. For x < — L, the pinned fluxon follows one of 
the two unstable manifolds of fixed point (arcsin7,0) of the reduced ODE (4). Similarly, for x > L 
the pinned fluxon follows one of the stable manifolds of the fixed point (arcsin7 + 27r, 0). Finally, for 
|x| < L the pinned fluxon corresponds to a part of one of the orbits of the phase portrait for the system 
with D — d. The freedom in the choice of the orbit in this system implies the existence of pinned 
fluxons for various lengths of the inhomogeneity. See Figure 1 for an example of the construction of a 
pinned fluxons when 7 = 0.15 and d — 0.2. Orbits of a Hamiltonian system can be characterised by 
the value of the Hamiltonian, hence there is a relation between the value of the Hamiltonian inside the 
inhomogeneity and the length of the inhomogeneity. The resulting pinned fluxon is in i7^(R) n C''^(M). 
As the ODE (4) usually implies that the second derivative of the pinned fluxon will be discontinuous, 
this is also the best possible function space for the pinned fluxon solutions. 

After analysing the existence of the pinned fluxons and having found a plethora of possible pinned 
fluxons when a bias current is applied to the Josephson junction (i.e., 7 7^ 0), we will consider their 
stability. First we will consider linear stability. To derive the linearised operator about a pinned 
fluxon 0pin(x; L, 7, d), write 0(x,t) = (/)pin(x; L, 7, d) + e^^i;(x, L, 7, d) and linearise about v — {) io 
get the eigenvalue problem 

>Cpin'^ = Ai;, where A = + «A, (6) 
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Figure 1: Phase portraits when 7 = 0.15 and d = 0.2. The dash-dotted red curves are the unstable 
manifolds of (arcsin7, 0), the dashed magenta curves are the stable manifolds of (27r + arcsin7, 0), and 
the solid blue curves are examples of orbits for the dynamcis inside the inhomogeneity. The bold green 
curve is an example of a pinned fluxon. 

and the linearisation operator Cpin{x] L, 7, d) is 



The natural domain for C^i^i is i/^(M). We call A an eigenvalue of >Cpin if there is a function v G i7^(R), 
which satisfies >Cpin(x; L, 7, d) = Av. This operator is self-adjoint, hence all eigenvalues will be real. 
Furthermore, it is a Sturm-Liouville operator, thus the Sobolev Embedding Theorem gives that the 
eigenfunctions are continuously differentiable functions in i7^(R). Sturm's Theorem [31] can be applied, 
leading to the fact that the eigenvalues are simple and bounded from above. Furthermore, if vi is an 
eigenfunction of Cpin with eigenvalue Ai and V2 is an eigenfunction of Cpin with eigenvalue A2 with 
Ai > A2, then there is at least one zero of V2 between any pair of zeros of vi (including the zeros at 
zboc). Hence, the eigenfunction vi has a fixed sign (no zeros) if and only if Ai is the largest eigenvalue 
of Cpin- The continuous spectrum of Cpin is determined by the system at ±00. A short calculation 
shows that the continuous spectrum is the interval (— oc, — y^l — 7^). 

If the largest eigenvalue A of Cpin is not positive or if Cpin does not have any eigenvalues, then 
the pinned fluxon is linearly stable, otherwise it is linearly unstable. This follows immediately from 
analysing the quadratic A = + aX. If A < 0, then both solutions A have non-positive real part. 
However, if A > is then there is a solution A with positive real part. Furthermore, the A- values of 
the continuous spectrum also have non-positive real part as the continuous spectrum of Cp[-a is on the 
negative real axis. 

The linear stability can be used to show nonlinear stability. The Josephson junction system without 
dissipation is Hamiltonian. Define P — (j)t^ u — {(j)^P)^ then the equations (1) can be written as a 
Hamiltonian dynamical system with dissipation on an infinite dimensional vector space of x-dependent 
functions, which is equivalent to H^{E}j x L^(R): 



/:pin(x;L,7,d) = D, 



XX 



L>cos(/)pin(x;L,7,d) = < 



Dxx ~ COS (/)pin (a;;L,7,d), \x\ > L; 
D^x - d cos(l)pin{x;L,-f,d), \x\ < L. 



(7) 
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and 



/oo |- 
P"^ + (t)l + 2 D{x; L, d) ( v/l - 7^ - cos (t>) dx 
-OO 



r»00 

7 ' 



roo rU 

/ [(/) — arcsin7 — 27t] + 7 / [0 — arcsin7] dx. 

Jo J -00 



(8) 



Here we have chosen the constants terms in the 7-integrals such that they are convergent for the fluxons. 
Furthermore, for any solution u{t) of (1), we have 

d 

—n{u) = P^dx < 0. (9) 

J —00 

As a pinned fluxon is a stationary solution, we have DT-L{(j)pi-^,0) = and the Hessian of T-L about a 
fluxon is 
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If Cpin has only strictly negative eigenvalues, then it follows immediately that ((/)pin, 0) is a minimum of 
the Hamiltonian and (9) gives that all solutions nearby the pinned fluxon will stay nearby the pinned 
fluxon, see also [10]. 

After this introduction, we will start the paper with an overview of simulations for the interaction 
of travelling fluxons and the inhomogeneity in (1) for various values of d, L, 7 and a. This will 
motivate the analysis of the existence and stability of the pinned fluxons in the following sections. 
We start the analysis of the existence and stability of pinned fluxons by looking at a microresistor 
with d = 0. The advantage of the case d = is that several explicit expressions can be derived 
and technical difficulties can be kept to a minimum, while it is also representative of the general case 
d < 1. It will be shown that for 7 = there is exactly one pinned fluxon for each length of the 
inhomogeneity. For 7 > 0, a plethora of solutions starts emerging. There is a minimum and maximum 
length outside which the inhomogeneity cannot sustain pinned fluxons. Between the minimal and the 
maximal length there are at least two pinned fluxons, often more. At each length between the minimum 
and maximum, there is exactly one stable pinned fluxon. If the length of the interval is (relatively) 
large, the stable pinned fluxons are non-monotonic. Note that stable non-monotonic fluxons are not 
possible in homogeneous systems, since for a homogeneous system the derivative of the fluxon is an 
eigenfunction for the eigenvalue zero of the operator associated with the linearisation about the fluxon. 
If the fluxon is non-monotonous, then this eigenfunction has zeros. As the linearisation operator is 
a Sturm-Liouville operator, this implies that the operator must have a positive eigenvalue as well, 
hence the non-monotonous fluxon is unstable. However, for inhomogeneous systems, the derivative 
of the fluxon is usually not differentiable, hence cannot give rise to an eigenvalue zero (since the 
eigenfunctions have to be C^) and stable non-monotonic fluxons are in principle possible. This shows 
that the inhomogeneity can give rise to qualitatively different fluxons. 

For the existence analysis of the pinned fluxons, the length of the inhomogeneity will be treated as 
a parameter. The pinned fluxons satisfy an inhomogeneous Hamiltonian ODE whose Hamiltonian is 
constant inside the inhomogeneity. It will be shown that the existence and type of pinned fluxons can 
be parametrised by the value of this Hamiltonian. The length of the inhomogeneity is determined by 
the value of the Hamiltonian and the type of pinned fluxon, leading to curves relating the length 2L 
and the value of the Hamiltonian inside the inhomogeneity. In [19] it is shown, in the general setting 
of an inhomogeneous wave equation, that changes in stability of the pinned fluxons can be associated 
with critical points of the length function relating L and the value of the Hamiltonian. The results of 
this paper together with Sturm-Liouville theory give the stability properties of the pinned fluxons in 
the general setting. 
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After giving full details for the case d = 0, for which the stability issue can be settled in dependent 
of [19], an overview of the results for d > is given. The general microresistor case (0 < d < 1) is 
very similar to the case d = 0. The microresonator case (d > 1) has some different features, but the 
same techniques as before can be used to analyse the existence and stability. We finish the analysis of 
the microresonator case by looking at the special case where microresonators approximate a localised 
inhomogeneity. We explicitly look at microresonators with d = ^ and L very small. For 7, a and /x 
small, the asymptotic results from [22] are recovered. Even in the limit of localised inhomogeneities, 
our work generalises [22], since our methods allows us to consider 7, a and ji larger as well. 

The paper concludes with some further observations, conclusions and ideas for future research. 

2 Simulations 

To put the analysis of the existence and stability of the pinned fluxons in the next sections in a wider 
context, we look first at simulations of the interaction of a travelling fluxon with an inhomogeneity. 
Recall that in absence of dissipation and induced currents (a = = 7), the system (1) without an 
inhomogeneity (i? = 1), has a family of travelling fluxon solutions (2), for each wave speed \v\ < 1, 
while if there is a small induced current and dissipation, but no inhomogeneity, then there is a unique 
travelling fluxon. 

First we look at the case o: = = 7 (no induced current, no dissipation) and the inhomogeneity 
of microresistor type with d = 0. If the length is too short, the fluxon will not be captured, but its 
speed will be reduced by the passage through the inhomogeneity. If the length of the inhomogeneity 
is sufficiently large, the travelling fluxon will be captured. Some radiation is released in this process 
and the fluxon "bounces" backwards and forwards around the defect, especially if the length is "just 
long enough". This is consistent with the results in [26] where a detailed analysis of the interaction 
of a fluxon with an inhomogeneity is studied in the case that no induced current and dissipation are 
present. An illustration is given in Figure 2. Note that the length of the defect which captures the 




Figure 2: Simulation of a travelling wave with speed v — 0.1 approaching an inhomogeneity with d = 
when there is no induced current (7 = 0) or dissipation (a = 0). The inhomogeneity is positioned in 
the middle (around the zero position). The length of the inhomogeneity on the left is 0.38 and the 
travelling fluxon is captured by the inhomogeneity; note that the "bounce" of the fluxon is a lot larger 
than the length of the inhomogeneity. The length of the inhomogeneity on the right is 0.36 and the 
pinned fluxon can just escape, but its speed is significantly reduced. 

fluxon is a lot smaller than the initial amplitude of the "bounce" of the fluxon. Observations suggest 
that the minimal length for the inhomogeneity to capture the travelling fluxon increases if the wave 
speed increases. 

Next we look at the system with a microresistor with d = 0, now with an induced current 7 = 0.1 
and varying lengths and values of a. We start again with an inhomogeneity of length 0.38 (L = 0.19). 
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When 7 = 0, this microresistor captures a fluxon with speed v = 0.1. With an induced current, it 
cannot capture a fluxon, however large we make a, i.e., however slow the fluxon becomes. This is 
illustrated in Figure 3. The microresistor slows the fluxon down for a while, but eventually the fluxon 




Figure 3: Simulation of a travelling fluxon approaching an inhomogeneity with d = when the induced 
current is 7 = 0.1. On the left, the length is 0.38. Here the dissipation is a = 0.9, but however large a 
is taken, the fluxon is never captured. In the middle and right plots, the length is 0.44. In the middle 
the dissipation is a = 0.48 and the fluxon is captured, on the right the dissipation is a = 0.47 and the 
fluxon can escape. 

escapes with the same speed as it had earlier (as this speed is unique in a system with a, 7 7^ 0). The 
simulations suggest that the smallest length which can capture a fluxon is 0.44 (L = 0.22). In the next 
section, it will be shown that for a, 7 7^ 0, there is a minimal length under which no pinned fluxon 
can exist. This explains why the inhomogeneity with the shortest length cannot capture even a very 
slow travelling fluxon. In Figure 3 it is illustrated that, if the length can sustain pinned fluxons, the 
capture depends on the dissipation (hence on the speed of the incoming fluxon). If the dissipation is 
sufficiently large, hence the speed sufficiently slow, the pinned fluxon will be captured. 

A longish defect in a microresistor will also capture the travelling wave and the resulting pinned 
fluxon is not monotonic, see Figure 4! The length of the inhomogeneity is substantial, so the stationary 




Figure 4: Simulation of a travelling fluxon approaching a longish inhomogeneity with d = when the 
induced current is 7 = 0.1 and dissipation is o: = 0.5. On the left, the length is 12.5, the travelling 
wave is captured and a non-monotonic pinned fluxon is formed. On the right, the length is 35 and the 
travelling wave escapes after a while, leaving in its wake a "bump" connecting 27v + arcsin7 at both 
ends. Note that the vertical scale and coloring is different in both figures; as a reference point, the 
travelling wave on the right is the same in both cases. 
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shape connecting the far field rest states at arcsin7 is a "bump". This "bump" is present at ah the 
rest states arcsin7 + 2k7v for 7 7^ as arcsin7 + 2k7v is not an equihbrium for the dynamics with 
d ^ 1. From a phase plane analysis it can be seen that the amplitude of the homoclinic connection to 
arcsin7 + 2k7v grows with the length L of the defect. As shown in Figure 4, for L = 6.25, the travelling 
fluxon travels into this "bump" and gets captured. The resulting pinned fluxon is not monotonic. In 
the next section, the family of all possible pinned fluxons is analysed and it is shown that for long 
lengths the stable pinned fluxon is non-monotonic. Moreover, it follows that there is an upper limit 
on the length of inhomogeneities that can sustain pinned fluxons. This is illustrated on the right in 
Figure 4. The travelling fluxon seems to be captured initially by the inhomogeneity, but after a while 
it escapes again. However large the dissipation is taken, this will always happen, illustrating that no 
pinned fluxons can exist. 

Next we consider a microresonator with d — 2. As before, we consider the case without an induced 
current (7 = 0) first. In this case, the fluxon is never captured. For the smaller lengths the fluxon 
reflects, for larger lengths the fluxon seems to get trapped, but it escapes after a while. This is 
illustrated in Figure 5 for a microresonator with length 0.1. In the next section, it will be will shown 
that a system with a microresonator and no induced current has indeed no stable pinned fluxons. 



( 




Figure 5: Simulation of a travelling wave approaching an inhomogeneity with d — 2 and length 0.1, 
when there is no induced current and no dissipation (7 = = a). The speed on the left is 0.21 and the 
travelling wave is bounced by the inhomogeneity. The length on the right is 0.22 and at first the pinned 
fluxon seems to be captured by the inhomogeneity, but after while it travels through the inhomogeneity 
and seems to resume its original speed. 

After the induction-less system, we consider a system with a microresonator with d — 2 and an 
induced current 7 = 0.1. As with the microresistor there is a minimum length, under which the 
microresonator cannot capture a fluxon. The simulations suggest that the minimum length is 0.42 
(L = 0.21). In Figure 6, it is illustrated that a microresonator with length 0.40 cannot capture a fluxon 
with Q{ = 0.9, whilst a microresonator with length 0.42 can capture a fluxon with a = 0.3, but it cannot 
for a = 0.29. This is consistent with the results in the next sections where it is shown that for a, 7 7^ 
there exists a minimal length under which no pinned fluxons can be sustained by the inhomogeneity. If 
the length can just sustain pinned fluxons, then there are both a stable and an unstable pinned fluxon 
close to each other. In the left panels of Figures 6 and 7 it can be observed that initially the travelling 
fluxon approaches the unstable pinned fluxon, but then reflects to the stable one and settles down. 

Finally we consider a microresonator with a longer length for which the travelling fluxon gets 
captured and becomes a non-monotonic pinned fluxon. In Figure 7, it is illustrated that for a microres- 
onator with d — 10, length 2 (L = 1), the travelling fluxon at 7 = 0.1 and a — 0.2 gets attracted to 
a non-monotonic pinned fluxon. Note that for microresonators (i.e., d > 1), the stable non-monotonic 
pinned fluxons have a "dip" as opposed to the ones for the microresistors which have a "bump" . 
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Figure 6: Simulation of a travelling wave approaching an inhomogeneity with d — 2 and length 0.1, 
when there an induced current (7 = 0.1). On the left and middle is a microresonator with length 0.42. 
On the left the dissipation is a = 0.3 and the fluxon is captured, whilst in the middle the dissipation 
is a = 0.29 and the fluxon escapes. On the right, the length is 0.4 and the dissipation is a = 0.9 and 
the fluxon still escapes as the length is too short for a pinned fluxon to exist. 
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Figure 7: Simulation of a travelling wave approaching an inhomogeneity with d = 10 and length 2, when 
the induced current is 7 = 0.1 and the dissipation is a — 0.2. The resulting wave is non-monotonic as 
can be seen on the right. Due to the weaker dissipation, it takes some time for the wave to converge to 
its stable shape. Initially, the travelling wave approaches the monotonic unstable pinned fluxon, then 
deflects from it and converges to the non-monotic stable one. 

3 No resistance (d=0) 

We now analyse the existence and stability of the pinned fluxons in a microresistor and a microresonator. 
First we consider the case when there is no resistance in the inhomogeneity, hence a microresistor with 
d — 0. This case provides a good illustration of the richness of the family of pinned fluxons, shows 
the essence of the analytic techniques for the existence and stability analysis, and has less technical 
complications than the more general values of d. The existence analysis for the case with no bias 
current (7 = 0) is quite different from the case when a bias current is applied (7 > 0). So we will 
consider them separately. 

3.1 Existence of pinned fluxons without applied bias current 

For 7 = 0, the pinned fluxon has to connect the stationary states at = and <j) — 2ti. In the 
background dynamics of the ODE (4) with = 1, the unstable manifold of (0,0) coincides with the 
stable manifold of (27r, 0), as follows immediately by analysing the Hamiltonian (5) with D = 1. These 
coinciding manifolds are denoted by a red curve in the phase portrait sketched in Figure 8. This curve 
represents the unperturbed sine-Gordon fluxon (2). The orbits generated by the Hamiltonian system 
with D = Q are straight lines. In Figure 8, samples of these orbits are given by the blue lines. Any 
blue line that crosses the red line can be used to form a pinned fluxon. An example is given in the 
panel on the right in Figure 8, where the green curve represents a pinned fluxon in n (7^(M). 

As can be seen from Figure 8, the value of the Hamiltonian inside the inhomogeneity is a convenient 
parameter to characterise the pinned fluxons. The points of intersection for the blue and red curves 
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Figure 8: Phase portraits of the ODE (4) for 7 = and d = 0. The red curve represents the coinciding 
stable and unstable manifolds of the asymptotic fixed points. The blue curves are orbits for the system 
inside the inhomogeneity. In the sketch on the right, the green curve represents a pinned fluxon. 

are denoted by ((/)in,Pin) respectively ((/>out5Pout) for the first respectively second intersection. It follows 
immediately that pi^ = :Pout and (j)o^t = 27r — (/)in. Furthermore, the expression for the Hamiltonian, (5), 
gives the following relations for (pi^ and pi^: = ^pfn ~ ~ ^os^in) {D = 1) and h = 
with < h < 2 where h is the value of the Hamiltonian inside the inhomogeneity. Thus 



Pm{h) — v2h and (l)in{h) — arccos(l — /i), with < h <2. 



Inside the inhomogeneity {\x\ < L), the pinned fluxon related to the value h satisfies h 
(j)x — \/2h. Hence the half-length L and the parameter h are related by 

L — I dx 



1^2 



J-L ^AnfW 4>X J<h:, 



TT — arccos(l — h) 



2h 



2h 



(10) 



thus 



(11) 



As the numerator is a monotonic decreasing function of h and the denominator is monotonic increasing, 
it follows immediately that L is a monotonic decreasing function of h. The function L takes values 
in [0, oc) as lim L{h) = oc and lim L{h) = 0. The h-L plot is given in Figure 9. We summarise the 
existence results for pinned fluxons without a bias current in the following lemma. 




Figure 9: Plot of the length L as a function of /i, the value of the Hamiltonian in the inhomogeneity, 
for 7 = and d = 0. 



Lemma 1 Let 7 = and d = 0. For any length 2L of the inhomogeneity, there is a unique pinned 
fluxon, for which the Hamiltonian inside the inhomogeneity has the value h(L) as implicitly given 
by (11). Define x* to he the shift such that (I)q{—L + x*) = (/)in (see (2) for the definition of (j)o), then 
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the pinned fluxon is given explicitly by 
(/)pin(x;L,0,0) = < 



0o(x + x*), 

TT— arccos(l— /i) 
I 



7r + 
00 (x 



X < — L, 
|x| < L, 
X > L. 



(12) 



3.2 Existence of pinned fluxons with bias current 

For 7 > 0, the pinned fluxon has to connect the stationary states at = arcsin7 and = 27r + arcsin7. 
In the background dynamics with D = 1 the unstable manifold of = arcsin7 coincides no longer with 
the stable manifold of 27v + arcsin7. Furthermore, the orbits of the dynamics inside the inhomogeneity 
are parabolic curves instead of straight lines. These two changes add substantial richness to the family 
of pinned fluxons. 

Let us first consider the phase portraits. In Figure 10 we consider 7 = 0.15 as a typical example 
to illustrate the ideas. In the dynamics with D = 1, the unstable manifolds to arcsin7 are denoted 
by red curves, while the stable manifolds to 27r + arcsin7 are denoted by magenta curves. The larger 
7 gets, the wider the gap between the unstable and stable manifold becomes. The dynamics within 





Figure 10: Phase portrait at 7 = 0.15 and d — 0. On the right is a zoom into the area around 
= (27r,0). 

the inhomogeneity with D = are denoted by blue curves. These blue curves are nested and can be 
parametrised with a parameter /i, using the Hamiltonian (5) with D = 0: 

where i?o(7) is given by the value of the Hamiltonian (5) on the magenta stable manifold {D = 1): 

^0(7) = 7arcsin7 — (1 — \/l — 7^) + 27r7. (13) 

Thus the value of h increases as the extremum of the blue curves is more to the right. 

For the existence of pinned fluxons, a blue curve has to connect the red unstable manifold with the 
magenta stable manifold. In Figure 10, the furthest left possible blue curve for which pinned fluxons 
may exist, is the one indicated with /i = 0. In the zoom on the right, it can be seen that this curve 
just touches the magenta stable manifold. The blue curve intersects the red unstable manifold twice, 
both points give rise to a pinned fluxon, as sketched in Figure 11. Obviously, the pinned fluxon in the 
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Figure 11: Phase portrait at 7 = 0.15 and d = with the furthest left blue curve for which pinned 
fluxons exist. There are two pinned fluxons possible, represented by green line. On the right is a zoom 
into the area around {(j)^(j)x) — (27r,0). 



second plot in Figure 11 will occur in a defect with a shorter length than the one in the first plot. The 
furthest right possible curve that gives rise to pinned fluxons is marked with /imax in Figure 10. This 
blue orbit touches the red unstable manifold and crosses the magenta unstable manifolds in 5 points. 
All these points represent different pinned fluxons, hence 5 pinned fluxons can be associated with this 
curve. Moreover, for h just below /imax, the blue curve intersects the red curve twice (while it still 
intersects the magenta curve 5 times: there are 10 different pinned fluxons associate to such value of h. 

In general, the pinned fluxons are determined by two points in the phase plane: the point where 
pinned fluxon enters the inhomogeneity (i.e. the crossing from the red unstable manifold to the blue 
orbit), this point will be denoted by ((/>in,Pin) and the point where the pinned fluxon leaves the inhomo- 
geneity (i.e. the crossing from the blue orbit to the magenta stable manifold), this point will be denoted 
by (^out^Pout)- Thus the points (0in,Pin) and (^out^Pout) are determined by the set of equations 

i7o(7) - 27r7 = " (1 " COS(/)in) + 70in, 



(14) 

2 ^out ~^ ^' 



i7o(7) = \pl^t - (1 - cos^out) + 70out. 
Combining the equations in (14), we get expressions for (/)in and (/)out: 

cos(/)in = 1 — (/i + 2ti^) and cos^out — h. (15) 
This is well-defined only if < /i < 2(1 — 7r7). Hence there are maximal values for 7 and /i, given by 

7max - and /i^ax = 2(1 - 7r7). 

TT 

If 7 > 7max5 then there is no blue curve that intersects both the red unstable orbit and the magenta 
stable orbit, hence no pinned fluxons exist if the applied bias current is larger than 7max- h> /imax, 
then the blue curve does not intersect the red manifold anymore. 

Furthermore, (/)in must lie on the red unstable manifold, hence arcsin7 < < 0max(7), where 
0max(7) is the maximal 0- value of the orbit homoclinic to arcsin7. As /i G [0, 2(1 — 7r7)], this implies 
that there are two possible values for and that > 0: 

(/)in = TT zb arccos(27r7 - (1 - h)) and pin = \/2 (i?o(7) + h - 7(/)in). 

Note that the unstable manifold left of arcsin7 only intersects with blue curves that have (^x < 0, hence 
those orbits can never connect to one of the stable manifolds of 27v + arcsin7. 
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The point (^out^Pout) has to lie on the magenta stable manifolds, so there can be up to five possible 
branches of solutions: 



1. 


V^out 


27r — arccosfl 


— }i \ with T>r^-,^^ ^ 

IV J \\ 11)11. />'out ^ ^5 


for all ^ h ^ h,^^^' 


2. 


0out 


= 27r + arccos(l 


- /i) with Pout > 0, 


for < /i < /i2 and pout < ioi h2 < h < h 


3. 


0out 


= 27r + arccos(l 


— h) with Pout ^ 


, for /l2 < < /imax; 


4. 


0out 


= 47r — arccos(l 


- h) with Pout > 0, 


for hi <h< /imax; 


5. 


0out 


= 47r — arccos(l 


- h) with Pout < 0, 


for hi < h < /imax. 



'max? 



Here /i2 is the /i-value such that the blue orbit intersects the magenta manifolds at the equilibrium 
(27r + arcsin7,0), i.e., /i2(7) = 1 — \/l — 7^, and /ii is such that the blue orbit touches the magenta 
manifold at (27r + 0max(7),O), the most-right point, thus /ii(7) = 1 — cos((/)max(7))- In all cases, 
bout I = x/2 (i7o(7) + /^-70out)- 

To satisfy /i2(7) ^ ^max(7)5 we need that 7 < 72 = i^^^T ^ 0.3104. If 7 > 72, then only pinned 
fluxons with (/)out = 27r ± arcsin7 and pout > exist. In order to have hi{^) < /imax(7)5 we need that 
7 < 71, where 71 is the implicit solution of cos (/)max(7i) + 1 = 27r7i, I.e., 71 

0.1811. If 7 > 71, then 

no pinned fluxons with (/)out = 47r — arcsin7 exist. On the intervals of common existence, we have 
< /i2(7) < ^1(7) < Vax(7), ^1(71) /imax(7i), ^2(72) /imax(72), scc Figure 12. 




0.4 



Figure 12: The extremal /i- values /ii(7), /?'2(7) and /imax(7)- 

In Figure 13, we have taken 7 = 0.15 and h = (/ii + /imax)/2 and have plotted all five possible 
pinned fluxons (i.e. all possibilities for (^out^Pout)) with (^[^ = tt — arccos(27r7 — (1 — h)). Obviously, 
five more pinned fluxons with the same (^out^Pout) are possible with (^[^ = tt + arccos(27r7 — (1 — h)). 

To determine the length of the inhomogeneity for the pinned fluxons, we use that on the orbits in the 
inhomogeneity (blue curves in the phase portrait) and (j)x are related by \(j)x\ = ^2 (i7o(7) -\- h — 70). 
Integrating this ODE, taking into account the sign of Pout, we get that the length of the pinned fluxons 
with Pout > is given by 



2L 



V2 



7 



\/H^ + h- 70in - \/Ho + h- 70out 
and for pout < 0, we have 

^/Ho^h- 7(/)in + Ho + h - 



7 



'out 



Pin - Pout 

7 



Pin - Pont 

7 



(16) 



(17) 
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Figure 13: The five pinned fluxons with (^[^ = tt — arccos(27r7 — (1 — h)) for 7 = 0.15, d = and 
h = (hi + /imax)/2. Note that only the pinned fluxon in the first panel is monotonic. In the L-h of 
Figure 14, the pinned fluxons in the first two panels are on the blue curve, the third one is on the red 
curve and the last two are on the green curve. 

These lengths are plotted in Figure 14 for 7 = 0.15. The blue curve is formed by the branches 1 and 
2, the red curve is branch 3 and the green curve is formed by the branches 4 and 5. This plot shows 
that there is a positive minimal and maximal length for the inhomogeneity to sustain pinned fluxons. 
Inhomogeneities with shorter or longer lengths will not be able to sustain pinned fluxons. Figure 14 
illustrates also that the maxima and minima of the possible length of the inhomogeneity are attained 
inside the interval (0, /imax), not at the endpoints. These extremal points will play an important role 
in the stability analysis as we will see in the next section. 

Remark 2 At 6 

out — 27r + arcsin7^ i.e h — h2, there is no bifurcation, although a solution ap- 
pears/disappears. To see that this is not a bifurcation point, we look at the disappearing solution in the 
limit h I For h \,h2, there is one solution with (/)out ^ 27r + arcsin7 and pout < and one solution 
with (/)out ^ 27r + arcsin7 and pont > 0. The solution with po^t < remains very close to 2ti + arcsin7 
for X > L. However, the solution with pout > is tracking almost all of the homoclinic connection to 
2ti + arcsin7. And in the limit h I h2 this solution ^'splits'^ into the pinned fluxon with pout — and a 
full homoclinic connection (fluxon- antifluxon pair). 

In general, the derivation of the existence of the pinned fluxons shows that for fixed 7 > and d = 0, 
there will always be a strictly positive minimal and maximal length for the existence of pinned fluxons. 
From Figure 12, it follows that the green curve of pinned fluxons with (/)out = 47r — arccos(l — h) is not 
present if 7 > 71. Similarly if 7 > 72, the red curve of pinned fluxons with (/)out = 27r + arccos(l — h) 
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Figure 14: The lengths of the pinned fluxons for 7 = 0.15 and d = 0. The lengths of the pinned 
fluxons with (/)out = 47r — arccos(l — h) are plotted in green (branches 4 and 5), the lengths of pinned 
fluxons with (/)out = 27r + arccos(l — h) and pout > are in red (branch 3). The lengths of the remaining 
pinned fluxons (branches 1 and 2) are indicated by the blue curves. The panels on the right zoom into 
the top and bottom and show that the minimal and maximal length are not obtained for /imax, but a 
smaller value. 

and Pout > are not present. Below we summarise the results for the existence of the pinned fluxons 
with an induced current: 

Theorem 3 For d = and every < ^ < ^, there are Lmin(7) ctnd Lmax(7); such that for every 
L G (i^min, ^max); there are at least two pinned fluxons (at least one for L =min or Lmaxy'- Furthermore 

limLmin(7) 0, limLmax(7) oc, 

7^0 7^0 

and 

lim Lmin(7) lim ^max(7) = a/ I ( arcsin ^ + Vtt^ - 1 J - W f ( arcsin ^ + Vtt^ - 1 - tt ) ^ 1.8. 

7tl/'7r Ttl/'TT V ^ ^ V ^ ^ 

For given L G [i^min, ^max]; the maximum possible number of simultaneously existing pinned fluxons 
is 6. For 7 > there exist no pinned fluxons. 

To relate the rich family of pinned fluxons which exists for 7 > with the unique pinned fluxons 
for 7 = 0, we have sketched the L-h curves for 7 = 0.001 in Figure 15. The bold blue curve converges 
clearly to the curve in Figure 9. These are the lengths associated with the pinned fluxons with (f)\^ — n — 
arccos(27r7— (1 — /i)) and (/)out — 27r— arccos(l — /i). Clearly there are some other convergent L-h curves as 
well. The length of the blue curve associated with the pinned fluxons with = 7r+arccos(27r7— (1 — /i)) 
and (j)out = 27r — arccos(l — h) goes to zero as expected. The convergent red and green curves can be 
associated with 47r-fluxons in the case 7 = 0. A 47r-fluxon is a connection between and 47r. This is 
not possible without an inhomogeneity, but with an inhomogeneity such connections are possible and 
some are stable. There are four possible 47r-fluxons if 7 = and the green and red curves converge to 
those solutions. For more details, see [18]. 
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Figure 15: L-h curves of the pinned fluxons for 7 = 0.001. 



3.3 Stability of the pinned fluxons with d = 

As seen in the introduction, the stabihty of the pinned fluxons is determined by the eigenvalues of the 
linearisation operator Cpi^i as defined in (7). For d = 0, the linearisation operator takes the form 



/:pin(x;L,7,0) 



Dxx - cos(/)pin(x;L,7,0), |x| > L; 



where 0pin is one of the pinned fluxons found in the previous section. 

When there is no induced current (7 = 0), expressions for the eigenvalues of Cpi^ can be found 
explicitly. Recall that for d = and 7 = 0, there is a unique pinned fluxon for each length L > 0, see 
Lemma 1. 

Lemma 4 For 7 = and d = 0, the linear operator >Cpin associated to the unique pinned fluxon in the 
defect with length L has a largest eigenvalue A^ax ^ (—1,0) given implicitly by the largest solution of 



-II 



V'2(l + COS 0in) + ^ (1 - cos 0in) = 

r 1 / \ (18) 

M+ W2(l + cos(/)in) tan /J"^'^ J- 

I ^ ^ \ \ y^2(l-cos(/)in) J 



where ji — \/l + A^ax ^ (0, 1) and the relation between (pi^v ct'^d L is given in (10) and (11). 

In Figure 16, A^ax is sketched as function of the half-length L of the pinned fluxon. The proof of 
Lemma 4 is quite technical; it is given in appendix A. 

Remark 5 For L large (hence 0in small), equation (18) has more solutions. Hence for those pinned 
fluxons >Cpin has some smaller eigenvalues in ( — 1,0) too. 

Corollary 6 // there is no induced bias current (^ — {)) and the microresistor has d — {), then the 
unique pinned fluxon in the defect with length L is linearly and nonlinearly stable. The pinned fluxon 
is asymptotically stable if a> 0. 

Next we consider the case that there is an induced bias current (7 > 0). In the previous section 
we have seen that in this case the pinned fluxons come in families, characterised by the blue, red and 
green curves in Figure 14. Locally along those curves, we can write either L as a function of h, or, 
/i as a function of L. Along those curves, we will look for changes of stability, i.e., find whether the 
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Figure 16: The largest eigenvalue of the linearised operator >Cpin at (i = and 7 = as function of the 
half-length L of the inhomogeneity. 

operator Cp[-a has an eigenvalue (recall that eigenvalues of Cp[n must be real. We will show that Cp[-a 
has an eigenvalue if and only if along the L-L curve we have ^ = or the pinned fluxon is isolated. 
Isolated pinned fluxons occur when 7 is maximal, i.e., 7 = ;^ or when 7 = 71, the maximal 7- value for 
which pinned fluxons with 0out = 47r — arccos(27r7 — 1) exist. This lemma is a special case of a more 
general theorem presented in [19]. The proof simplifies considerably in this case. 

Lemma 7 For any 7 > 0^ the linear operator Cpin(x; L, 7, 0) has an eigenvalue zero if and only if 

• dh 

• or ^ — ^ (this eigenvalue zero is the largest eigenvalue); 

• or 7 = 71 ^ 0.18; the solution of 

cos 0niax(Ti) ~^ 1 — 27r7i (see section 3.2), and (f)p\n is such that 
0in = 0out = 47r — arccos(27r7i — 1) = 27r + 0max(7i) (this eigenvalue zero is not the largest 
eigenvalue). 

Proof First we observe that differentiating (4) with respect to x shows that (j)'^:^^ satisfies ^pin^pin — 0- 
However, it follows immediately from (4) that ^p-^^ is not continuously different iable, except when there 
exist A;± G N such that (j)^^ — k-7v and (^out — k^Tv. From the existence results, it follows that this 
happens only if 7 = In this case, there is only one pinned fluxon and the blue curve in Figure 14 
has become a single point (there are no red or green curves). 

In all other cases, 0^-^^ C-'^(R) D i7^(R) so 0^-^^ is not an eigenfunction with the eigenvalue 
zero. However, 0^-^^ still plays a role in the eigenfunction related to any eigenvalue zero. Indeed, 
on both intervals (oc, — L) and (L,oc), the second order linear ODE >Cpin'0 = has two linearly 
independent solutions. As the asymptotic system is hyperbolic, one solution is exponentially decaying 
whilst the other is exponentially growing. Thus if the linear operator C has an eigenvalue zero, then 
the eigenfunction in the intervals (— oc, L) and (L, oc) must be a multiple of the exponentially decaying 
solution. As ^p-^^ is exponentially decaying for |x| ^ oc and satisfies >C0p-j^ = for |x| > L, it follows 
that for any eigenvalue zero, the eigenfunction must be a multiple of c/)^^^^ for |x| > L, unless c/)^^^^ = 0. 
The case c/)^^^^ = happens only when (/)out = 27r -h arcsin7 and x > L. In this case, the appropriate 

eigenfunction is a multiple of e~ 
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Next we look inside the inhomogeneity, i.e., |x| < L. The hnearised problem inside the defect for 
an eigenvalue zero can be solved explicitly and gives an eigenfunction of the form A + B{x + L)^ with 
A and B free parameters and |x| < L. 

To conclude, if the linear operator C^i^ has an eigenvalue zero, and (/)out ^ 27r + arcsin7 (we will 
consider this case later), then the eigenfunction is of the form 



A + S(x + L), |x| < L, 



where A, B and K are free parameters. We have to choose the free parameters such that '0 is con- 
tinuously different iable at zbL. As there are only three free parameters and four matching conditions, 
this will give us a selection criterion on the length L for which an eigenvalue zero exists. The matching 
conditions are 

A^^^iJ-L-), B^<j,';,^{-L-), B^K<j,';^{L^), and A^2BL ^ K<j,'^,^{L^), 

where the notation (l)'^;^{-L-) = lim^^_L 0pin(^), 0pin(^^) = lim^^L 0pin(^), etc. Using that Pin/out = 
0p-j^(^L) and 7 + (/)p-j^(zbL^) = sin(/)(±L) = sin ^j^/out 5 ^his can be written as 

A = Pin, B = sin(/)in - 7, B ^ K(sin(/)out - t), and A + 2SL = Kpout- 
Equations (16) and (17) show that L — hence the parameters are given by 

A = Pin, B = sin 0in - 7, and K(sin (j)out - 7) = sin (/)in - 7 
and the compatibility condition on L, or equivalently /i, is 

= Pin sin0in(sin(/)out - 7) - Pout sin 0out (siu (/)in - 7). (19) 

To derive this expression, we have multiplied the remaining equation [A+2BL = A^Pout] with 7(sin(/)out — 
7). This term would be zero if sin^out — 7? hence (/)out = 27r + arcsin7 but this case is not considered 
now. 

For completeness, we also consider the case where we assume that the eigenfunction vanishes for 
X < —L. If this is the case, then matching at x = —L gives immediately that A = = B. Thus this 
leads to a non-trivial eigenfunction only if (/)p-j^(L) = = lim^^^i, (/)p-j^(x). In other words, when (/)pin is a 
fixed point for x > L. This happens only if 0out = 27r + arcsin7. This case we will be considered later. 

Next we link the expression (19) to the derivative of L with respect to h. As L = P^n-Pout ^ ^j^^ 
derivatives of pin and pout are needed. Differentiating (14) and (15), we get 

Pin = 1 - 70in(^)^ sm^in = 1 and pout = 1 - 70out(^)^ sm^out = 1- 

Thus differentiating L — ^'''27°''^ gives that 

Pin sin (/)in Pout sin (j)ont 37" = ^ [Pout sin 0out (sin (/)in - 7) - Pin sin (/)in(sin (/)out - 7)] (20) 
ari Z7 

So we have shown that if 0out 7^ 27r + arcsin7 and the operator >Cpin has an eigenvalue zero, then either 
^(/i,7) = or Pin sin (/)in Pout sin 0out = 0. Considering pin sin 0in Pout sin 0out = in more detail, we 
get: 

• sin^out — would mean that (/)out — 27r. Going back to the compatibility condition (19), this 
implies that 7PinSin(/)in = 0, which only happens if also sin^in = or pin = 0. In the existence 
section we have seen pin > 0, hence 7Pin sin^in = can only happen if — tt, hence if 7 = -; 
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• sin(/)in = implies that (^[-a — tv. Going back to the compatibihty condition (19), this imphes 
that 7Pout sin^out — 0, which only happens if also sin^out = or pout — 0. Hence either "y — ^ 
or 7 = 7i, as the case (^out = 27r + arcsin7 is excluded at this moment; 

• l^in 7^ as we have seen before; 

• Pout — happens if (/)out = 27r + arcsin7 or (/)out = 27r + 0max(7)- Going back to the compatibility 
condition (19), this implies that sin(/)in(sin0out ~ t) = 0. Since tt — arcsin7 < 0max(7) < 27r, 
this implies this only happens if sin^in = 0, which case is considered before. 

So altogether we have if (^out ^ 27r + arcsin7 and the operator Cpi^i has an eigenvalue zero, then either 

• 0in = and 0out — 27r, which only happens when 7 = The eigenfunction in this case is ^pj^' 
which does not have any zeros, hence the eigenvalue zero is the largest eigenvalue. 

• 0in ~ ^ and (t>out = 27r + 0max(7) (i-e. Pout = 0), which only happens if 7 = 71. In this case the 

eierenfunction is 6'- for x < L and — ? — 7 6^- for x > L. This eierenfunction has a zero at 

X = L, hence the eigenvalue zero is not the largest eigenvalue. Note that when 7 = 71 the green 
L(h) curve in Figure 14 has degenerated to an isolated point related to the pinned fluxon 
considered in this case. 

To show that the converse is true, we look at the three cases f^(^5 7) = 0, 7 = and 7 = 71 and 
(^out^Pout) = (27r + 0max, 0). It is Straightforward to verify that the eigenfunctions as described earlier 
can be constructed in those cases. 



Finally we look at the case 0out = 27r + arcsin7. In this case 7 < and h — h2 — 1 — \/l — 7^. 

Furthermore, the pinned fluxons satisfies (j)'^:^^ = for x > L. In this case, the general form of an 
eigenfunction for an eigenvalue zero is 



A + B {x -\- L), \x\ < L, 



where A, B and K are free parameters. We have to choose the free parameters such that is contin- 
uously differentiable at = ±L, i.e. 



A^^^iJ-L-), S = (/);;iJ-L-), K = A + 2BL, and B^-K^/l^. 

As L = P^n-Pout ^ ^j^jg implies that A = pin, S = sin^in - 7 and K = pir^^^^ ^ith the matching 
condition 



7(sin(/)in -7) = - v/l - 7^ sin(/)inPin. (21) 

If 0in = + arccos(27r7 — y^l — 7^), then sin^in < and (21) cannot be satisfied as pi^ > and 
7 > 0. If (/)in = TT — arccos(27r7 — \/l — 7^), then the phase portrait in the existence section shows 
that sin(/)in > 7, thus (sin^in — 7) > and again (21) cannot be satisfied. Thus no eigenvalue zero can 
occur at (pout = 27r + arcsin7. □ 

Lemma 7 allows us to conclude the stability of pinned fluxons. 

Theorem 8 For d = 0^ every < j < and every L G [i^min(7), ^max(7)]; there is exactly one 
stable pinned fluxon. This pinned fluxon is linearly and nonlinearly stable (and asymptotically stable 



for a > 0). For L sufficiently large > y 7r+arcsin7+arccQs(27r7 y^i 7 ^ j ^ stable pinned fluxons are 
non-monotonic. 
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Figure 17: Stability for (i = and 7 = 0.15. The bold magenta curve represents stable solutions, all 
other solutions are unstable. On the right there is an example of a stable monotonic pinned fluxon (at 
L — 0.38) and a stable non-monotonic one (at L — 10), Both stable pinned fluxons have /i = 1, i.e. 
they are near minimal respectively maximal length, which are at Lmin — 0.35 and -t/max — 

10.13. 

See Figure 17 for an illustration of this theorem. 



Proof If 7 = then only the inhomogeneity with half-length exactly L = W | ( arcsin ^ + a/tt^ — 1 j — 



arcsin - + a/tt^ — 1 — vr) p:^ 1.8 has a pinned fluxon. From Lemma 7, it follows that the linearisa- 



tion for this pinned fluxon has a largest eigenvalue 0, so this pinned fluxon is linearly stable. 

In Corollary 6, we have seen that the unique pinned fluxons for 7 = are stable. 

If < 7 < then there are at least two pinned fluxons if L E (i^min, ^max), see Theorem 3. As 
seen before, the L-h curves for the pinned fluxons form three isolated curves: 0out = 47r — arccos(l — h) 
(green curve), the (red) curve of pinned fluxons with 0out = 27r + arccos(l — h) and pout > (exists 
for h > /i2), and the other pinned fluxons (blue curve). The colour coding refers to Figures 14 and 17. 
The fluxons on the blue curve exist for all < 7 < ^; the existence of the other curves depends on the 
value of 7. 

The linearisation about the pinned fluxon at the minimum on the red curve has an eigenvalue zero. 
At this point, the associated eigenfunction is a multiple of 0^-^^ for x > L. On the red curve, pout > 
and 0pin(^) < for X large. Thus this eigenfunction has a zero. Using Sturm-Liouville theory, we can 
conclude that the eigenvalue zero is not the largest eigenvalue. As there is only one fluxon with = 
on the red curves, all pinned fluxons on the red curve are linearly unstable. 

Similarly, the minimum and maximum on the green curve are associated with pinned fluxons whose 
linearisation has an eigenvalue zero. Again, the associated eigenfunction for x > L is a multiple of 
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0p-j^. As for the red curve, at the minimum we have pout > and 0pin(^) < for x large. Thus this 
eigenfunction has a zero and we can conclude that the eigenvalue zero is not the largest eigenvalue. 
The green curve is a closed curve with only two points with ^ = 0, so the eigenvalue zero at the 
maximum cannot be the largest eigenvalue either. So we can conclude that all pinned fluxons on the 
green curve are linearly unstable. 

Finally we consider the blue curve. We use the stability of the pinned fluxons at d = 0, 7 = to get 
a conclusion about the stability of the pinned fluxons on this curve. The solutions that can be continued 
to 7 = are the connections between — n — arccos(27r7 — 1 + /i) and (/)out = 27r — arccos(l — h). 
Hence those solutions are stable. Now using that zero eigenvalues can only occur if L{h) has a critical 
point, the blue curve can be divided in stable and unstable solutions. The stable solutions are the part 
of the curve L{h) curve between the minimum and maximum that contains the pinned fluxons with 
0in = TT — arccos(27r7 — 1 + /i) and (/)out = 27r — arccos(l — h). The pinned fluxons in the other part are 
unstable as the zero eigenvalue is simple. It can be verifled that the eigenfunctions related to the zero 
eigenvalues on this curve do not have any zeroes indeed. 

So altogether we can conclude that for each length there is exactly one stable and at least one 
unstable solution. The stable fluxons are non-monotonic if L is larger than the length of the fluxon 
at /i = /i2(7) = 1 — a/I — 7^ with (j)^^ — n — arccos(27r7 — y^l — 7^) and (/)out = 27r + arcsin7, hence 

^ / 7r+arcsin7+arccos(27r7— y^l^^T^) l— . 

4 General case {d > 0) 

After analysing the existence and stability of pinned fluxons in microresisors with the d = in full 
detail, in this section we will sketch the existence and stability of the pinned fluxons for a general 
microresistor or microresonator. 



4.1 Microresistors (0 < < 1) 

The existence of pinned fluxons for < < 1 follows from similar arguments as for the case = 0. 
Using the matching of appropriate solutions in the phase planes again, it can be shown that pinned 
fluxons exist for < 7 < The Hamiltonian dynamics in the inhomogeneity satisfles the relation 

- d{l - cos 0) + 70 = ^0(7) + h, 

where /i is a parameter for the value of the Hamiltonian as before. The case 7 = (no induced current) 
is more or less identical to before, with a unique pinned fluxon for any L > 0. For 7 > 0, a similar 
calculation as in the case d = shows that there are two possible entry angles: 

(Pin — 7T — arccos I — I — 7T + arccos I — I 

and up to three possible exit angles: 

(/)out = 27r - arccos (^^z^^) , 0out = 27r + arccos (^^z^^) , or (/)out = 47r - arccos (^r^) ' 

with < h < 2(1 — d — 7r7). If 7 > d > (i.e., d is sufliciently close to zero), then there is still a minimal 
length ivmin(7) > and a maximal length Lmax(7) foi" the inhomogeneity at which pinned fluxons can 
exist. However, if 7 is less than d (0 < 7 < d), then there is no upper bound on the possible length of 
the inhomogeneity anymore, i.e., Lmax = oc. This new phenomenon appears for 7/d < 1, due to the 
fact that now the dynamics in the inhomogeneity have flxed points at (0,p) = {2k7v + arcsin(7/d), 0), 
A; G Z. If /i corresponds to an orbit which contains such a flxed point, then the length of an orbit 
with Pout < goes to inflnity. To illustrate this, in Figure 18, we have sketched the phase portraits for 
d = 0.2 and 7 = 0.15 < d and 7 = 0.22 > d. 
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Figure 18: Phase portrait at d = 0.2 and 7 = 0.15 (left) and 7 = 0.22 (right). Note that in the left 
graph, the third blue orbit has a fixed point. So the pinned fluxon with cpo^t = 27r + arccos ^ "'^Y-d^ ) 
and Pout < does not exist for this /i- value. Nearby pinned fluxons will be in a defect with a length 
that goes to infinity. In the right graph, there are no fixed points anymore as ^ > 1. Thus the defect 
lengths for which pinned fluxons exist are bounded. 



As before, the length of the inhomogeneity for the pinned fluxons parametrised with h can be 
determined by using the relation \(j)x\ — \/2(i7o(7) h-\- d{l — cose/)) — 7^) and integrating the ODE, 
taking care of the sign of (f)x- The resulting integrals cannot be expressed analytically in elementary 
functions anymore, but they can be evaluated numerically. To illustrate this, we have determined the 
L-h curves as function of /i for d = 0.2 and 7 = 0.15 (7 < d) and 7 = 0.22 (7 > d). The L-h curves are 
presented in Figure 19. Note the unbounded length curve for 7 = 0.15. 
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Figure 19: L-h curves for d = 0.2 and 7 = 0.15 (left) and 7 = 0.22 (right). For 7 = 0.15, the L-h 
curves are unbounded as 7 is less than d. The color coding is as before, hence the bold magenta curve 
correspond to the stable fluxons. 

In the following theorem, we summarise the existence of pinned fluxons for < d < 1 and give their 
stability. 

Theorem 9 For < d < 1 and 

• 7 = 0^ there is a unique stable pinned fluxon for each L > 0; 
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• < 7 ^ min (d, , there is a minimal length Lmin(7) > such that for all L > Lmin there 
exists at least two pinned fluxons (one for L = i^min)- For each L > Lmm; there is exactly one 
stable pinned fluxon; 

• d < 7 < there are minimal and maximal lengths, ^min(7) > respectively Lmax(7) such 
that for all Lmin < L < L^ax there exists at least two pinned fluxons, one pinned fluxon if L is 
maximal or minimal, and no pinned fluxons exist for other lengths. For each Lmin ^ L < L^ax; 
there is exactly one stable pinned fluxon; 

• for 7 > there exist no pinned fluxons. 

Note that the third case wih be relevant only if < d < 

To prove the stability result for the pinned fluxons, we will use Theorem 3.1 from [19]. In [19], 
the stability of fronts or solitary waves in a wave equation with an inhomogeneous nonlinearity is 
considered. It links the existence of an eigenvalue zero of the linearisation with critical points of the 
L-h curve. The proof has similarities with the proof of the case d = in Lemma 7, but several extra 
issues have to be overcome. Theorem 3.1 of [19], applied to our pinned fluxons for < d < 1, leads to 
the following lemma, which is very similar to Lemma 7 which holds for the microresistor with d = 0. 

Lemma 10 IfO<d<l, then the linear operator >Cpin(x; L, 7, d) has an eigenvalue zero if and only if 

• dh 

• or ^ — (this eigenvalue zero is the largest eigenvalue); 

• or ^ is such that it solves (1 — d)(cos 0max(7) + 1) = 27r7 and the pinned fluxon is such that 
0in = 27r + 0max(7) (this eigenvalue zero is not the largest eigenvalue). 

The verification of Lemma 10 can be found in [19, §3.5]. As far as the special cases in this lemma is 
concerned, if 7 = or 7 is such that it solves (1 — d)(cos(0max(7) + 1) = 27r7 and the pinned fluxon 
is such that (/)in = 27r + 0max(7), then the pinned fluxon under consideration corresponds an isolated 
"green" point and ^ does not exist. In the case of 7 = there is exactly one value of the length L 
for which there exists a pinned fluxon. In the other case, there are more pinned fluxons, but on other 
branches. In the case of an isolated pinned fluxon, either the derivative of the pinned fluxon is an 
eigenfunction with the eigenvalue zero or a combination of multiples of the derivative of the pinned 
fluxon is an eigenfunction. 

The stability result of Theorem 9 follows by combining Lemmas 4 and 10. 

Proof of Theorem 9 The existence is described at the first part of this section, in this proof we 
focus on the stability. For < d < 1 and 7 = 0, there is a unique pinned fluxon for each length L. It 
is straightforward to show that for each < d < 1, the length function L{h) is monotonic decreasing 
m h. Thus ^ ^ and none of the pinned fluxons has an eigenvalue zero. As all pinned fluxons are 
nonlinearly stable for d = (Lemma 4) and no change of stability can happen, all pinned fluxons with 
7 = are nonlinearly stable for all < d < 1. 

If < d < 1 and < 7 < ^— ^, then the L-h curve follows as a smooth deformation from the curve 
for d = 0. And the unique stable pinned fluxon for each length follows. 

If < d < 1 and 7 = then the pinned fluxon is an isolated point and Lemma 10 gives that it 
is stable. □ 
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4.2 Microresonator {d > 1) 

The existence results of pinned fluxons for d > 1 are slightly different from the ones for d < 1. The main 
difference is the type of solutions used in the inhomogeneous system. For d < 1, we used solutions that 
were part of unbounded orbits or homo/heteroclinic orbits in the phase plane. For d > 1, we have to use 
periodic orbits. The most simple way to understand this crucial difference between the microresistor 
and the microresonator case is to consider the phase portraits without applied bias current (7 = 0)- 
see Figure 20. When d < 1, respectively d > 1, the (red) heteroclinic orbit of the system outside the 




Figure 20: Phase portraits at 7 = for various values of d. The red curve is the heteroclinic connection 
at d = 1. The blue curves are orbits for d — ^ and the green ones are orbits for d — 2. 

inhomogeneity is outside, resp. inside, the (blue resp. green) heteroclinic orbit of the system inside 
the inhomogeneity - see Figure 20. As a consequence, a pinned defect can only be constructed with 
(unbounded) orbits that are outside the (blue) inhomogeneous heteroclinic orbit in the microresistor 
case, while one has to use bounded, periodic orbits in microresonator case - see the green lines in 
Figure 20. 

One consequence is that if one solution for a inhomogeneity of a certain length exists, then there 
are also solutions for inhomogeneities with lengths that are this length plus a multiple of the length of 
the periodic orbit. This implies that the number of pinned fluxons for a defect of length L may grow 
without bound as L increases - which is very different from the microresistor (d < 1). We will focus on 
the existence of solutions which use less than a full periodic orbit as the other ones follow immediately 
from this. 

Using similar techniques as in the previous sections, it can be shown that if d is the solution of 
— ^ + arcsin ^ + \/d^ — 1 + d— 1 = 0, (d 4.37), then for d > d, pinned fluxons exist for any < 7 < 1. 
If d < d, then pinned fluxons exist for < 7 < 7max5 where 7max(rf) is the (implicit) solution of 
— 27r7 — 7 (arcsin 7 — arcsin 2) + A/d^ — 7^ — y^l — 7^ + (d — 1) = 0. 

For illustration, phase portraits for d = 4 and various values of 7 are sketched in Figure 21. This 
illustrates that the solutions used in the inhomogeneous system (blue lines) are all part of a periodic 
orbit. Note that for 7 > both unstable manifolds of arcsin 7 and only the unbounded stable manifold 
of 27r + arcsin 7 are used as opposed to the microresistor case where only the bounded unstable manifold 
of arcsin 7 and both stable manifolds of 2ti + arcsin 7 are used. 

As before, the dynamics in the inhomogeneity satisfies the relation 

- d(l - cos 0) + 70 = i7o(7) + 
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Figure 21: Phase portrait at d = 4 and 7 = 0.2 (left), 7 = 0.5 (middle) and 7 = 0.95 (right). As before, 
the red curves are the unstable manifolds to arcsin7 and the magenta ones are the stable manifolds to 
27r + arcsin7. The blue curves are orbits inside the inhomogeneity. The inner blue curve with angles 
between and 27r is the curve with the minimal /i- value for which pinned fluxons exist. The blue curves 
can continue to relate to pinned fluxons up to (but not including) the blue homoclinic connection to 
arcsin(^). Some of the periodic orbits with negative angles will also play a role in the construction of 
the pinned fluxons. If 7 = 0.95 > 7max(rf) ^ 0.9 (right plot), the blue homoclinic orbit (that encloses 
the red limit state) does not intersect the magenta stable manifold; illustrating that there cannot be 
pinned fluxons for 7 > 7max(<i)- 

where /i is a parameter for the value of the Hamiltonian. Again it can be shown that the entry and 
exit angles satisfy 



where now —2{d — 1) < h < /imax- Here /imax corresponds to the /i- value of the orbit homoclinic to 
arcsin ^ in the inhomogeneous system; it can be shown that /imax < 0. As we use periodic orbits 
inside the inhomogeneity, the entry and exit angles will differ by less than 27r. For any h value in 
[-2(d-l),/i 

max)? there will be pinned fluxons with entry angles between arcsin ^ and 27r. For 7 small 
relative to d, entry angles less than arcsin ^ are also possible and they can be related to smaller (more 
negative) h values. The value for the exit points is always positive, while the entry points can have 
both positive and negative values if the entry angle is larger than arcsin^. The pinned fluxons 
with entry angles less than arcsin 7 have only negative pin- values and hence those pinned fluxons are 
non-monotonic and "dip down" . 

For 7 = 0, at least one pinned fluxon exists for each L > 0. If L is sufficiently large, there will 
be more pinned fluxons. This is different to the case with d < 1, where for 7 = 0, there is a unique 
pinned fluxon for each length, it is due to the fact that the pinned fluxons are buiklt from periodic 
orbits (that may be travelled in various waus before leaving the inhomogeneity). For 7 > 0, there is 
minimum length Lmin such that there are at least two pinned fluxons for each length L > Lmin (one 
for L minimal). The L-h curves for d — A and various 7 values are given in Figure 22. Only lengths of 
the pinned fluxons that use less than a full periodic orbit are plotted. 

In the following theorem, we summarise the existence of pinned fluxons for d > 1 and give their 
stability. 

Theorem 11 Let d he the solution of — ^ + arcsin ^ + \/d'^ — 1 + 1 = (d ^ 4.37) and for d > 1, let 
7max(rf) be the (implicit) solution of — 27r7 — 7 (arcsin 7 — arcsin ^) + — 7^ — \/l — 7^ + (d — 1) = 0. 

• For d > 1 and 7 = 0; there is at least one pinned fluxon for each L > and all pinned fluxons 
are unstable; 

• For 1 < d < d and < 7 < 7max(rf); there is a minimal length Lmin(7) > such that for all 



cos (/)[ 



27r7 + d-l + h 
d-1 
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Figure 22: L-h curves at d = 4 and 7 = (left), 7 = 0.2 (middle) and 7 = 0.5 (right). The blue and 
red curves are associated with pinned fluxons with arcsin7/(i < (^[^ < tt. The red curve are pinned 
fluxons with < and arcsin7 < (^[^ < tt. The pinned fluxons in the blue curve have pi^ > for 
(/)in > arcsin7 and j)in < for (/)in < arcsin7. The green curves are associated with pinned fluxons with 
0in > TT. In the middle panel (7 = 0.2) there are also black curves, which are associated with pinned 
fluxons with (f)[^ < arcsin ^. The solid black curves are lengths for pinned fluxons with — 27r + arcsin ^ < 
0in < arcsin ^, the dashed ones for pinned fluxons with — 47r + arcsin ^ < < —2ti + arcsin 2, the 
dotted ones for pinned fluxons — Gtt + arcsin ^ < (f)\^ < — 47r + arcsin ^. 

L > Lmin there exist at least two pinned fluxons (one for L = i^min)- For each L > Lmm; there is 
at least one stable pinned fluxon. 

• For d > d and < 7 < 1^ there is a minimal length Lmin(7) > such that for all L > Lmm there 
exist at least two pinned fluxons (one for L = i^min)- For each L > Lmm; there is at least one 
stable pinned fluxon. 

In Figure 22, the stable pinned fluxons are the pinned fluxons on the increasing part of the lower 
right blue curve. Note that these pinned fluxons are non-monotonic past the meeting point with the 
red curve, hence for most lengths. The fluxons on the other blue curve and red and green curves 
are unstable. As before, the proof of the stability properties of Theorem 11 is based on Theorem 3.1 
from [19]. The proof of Theorem 11 is very similar to the proof of Theorems 8 and 9. The main 
difference is that we can not track our stability arguments back to the case d = (i.e. Lemma 4) as 
we did before. The role of Lemma 4 will now be taken over by Lemma 12 in Appendix A, in which it 
is explicitly established tht the pinned fluxon on the blue curve has exactly one positive eigenvalue for 
7 = and d near one. 

The stability of the fluxons on the black curves can not easily be related to fluxons at 7 = (they 
"split" in a homochnic "dip" and a fluxon for 7 = 0). So a stability analysis for this case goes outside 
the scope of this paper. In section 5, we will show numerically that there are some stable fluxons on 
the black curve. 
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4.3 A microresonator approximating a localised inhomogeneity 



There have been quite a number of investigations on the influence of a locahsed inhomogeneity, i.e., 
D{x) = (1 + iJi5{x)) or D{x) = (1 + Yl^=i l^i^i^ ~ ^i)) (!)• this section we wifl conflrm that our 
existence and stabihty results, apphed to short microresonators with large d, reproduce in the limit 
for L ^ and d ^ oo the existence and stability results for pinning by microshorts in [22]. In [22] it 
is shown that for D{x) = (1 + iJi5{x)) and 7, /x, and a of order with e small and ^ < + 0{e)^ 
there are one stable and one unstable pinned fluxon, both approximated by 0o(^ — -^0) + where 
Xo are the two solutions of — ^ + sech^XtanhX. 

To approximate the localised inhomogeneities of 5-function type with flnite length ones, we look at 
microresonators with length L — 1/(2^) and d = 1 + /xt) for D large. Thus the microresonators have 
short lengths and we can restrict to pinned fluxons with 

/2^7 + d-l + /i\ / , 2^7 + /i\ 
0in = arccos — ; = arccos 1 H , Pin > 0, 

and 

fd-l + h\ / h 

0out — arccos — — arccos 1 H 

\ d — I J \ jjA) 

Hence the pinned fluxons of [22] correspond to solutions on the lower blue curve in Figure 22. Intro- 
ducing /i = /xc)[), we get that — 2 < [} < and we are interested in [} away from as [)- values close to 
zero correspond to long lengths. Using the expressions for (/)in, Pin ^out? and pout and the ODE for the 
pinned fluxon, we can derive an asymptotic expression for the length L([}) if ^ is large and 7, ji are 
order where e is small: 



Thus L(f)) has a minimum at f) = — | + + e) and the condition L(f)) = l/2i) can be satisfied 

if the cubic \i^{2 + f)) = 2^ + 0{<3-^ + e) can be solved for some f) < 0. For f) < 0, this cubic has 

fi 

a maximum at f) = — | + 0{l)~^ + £), thus L([)) = ^ has two solutions with [} between —2 and iff 
^ < + (9(c)~-'^ +s) (i.e., there are no solutions for 7//X too large). From the analysis in the previous 
section, we can conclude that this corresponds to one stable pinned fluxon (least negative value of [)) 
and one unstable pinned fluxon. 

Finally, for 7 = 0{e)^ with e small, both the unstable manifold to arcsin7 and the stable manifold 
to 2ti + arcsin7 are close to the heteroclinic connection for the unperturbed sine-Gordon equation. 
Thus for X > L, we have — (j)o{x — Xq) + 0{e)^ where 00 is the shape of the stationary fluxon 

in the sine-Gordon equation (and a similar relation for x < —L). Substituting this into the equation 
for 0out5 with L = 1/2^ (hence [) is a solution of the cubic introduced earlier), we get that Xq is one of 
the two solutions of — ^ + sech^XtanhX = 0. 



5 Conclusions and further work 

This paper exhibits a full analysis for the existence and stability of pinned fluxons in microresistors 
and microresonators for which the Josephson tunneling critical current is modelled by a step-function. 
It is shown that for flxed d (Josephson tunneling critical current inside the inhomogeneity) and flxed 
bias current 7 > 0, there is an interval of lengths for which a rich family of pinned fluxons exists. In the 
case when an induced current is present, there is a lower bound on the length of inhomogeneities for 
which pinned fluxons can exist. If the inhomogeneity is too short, no pinned fluxons can be sustained. 
The lower bound on the length increases if the induced current increases. For microresistors with a 
sufficiently large induced current, there is also an upper bound on the length for pinned ffuxons and 
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the upper and lower bounds collide when the maximal value of the induced current for which pinned 
fluxons can exist, is attained 

Compared to the case of homogeneous wave equations, a new phenomenon is observed: longer 
microresistors and microresonators have non-monotonic stable pinned fluxons. In the case of mi- 
croresistors {d < 1), the non-monotonic stable pinned fluxons have a "bump" inside and behind the 
inhomogeneity and the values in the bump exceed the asymptotic state 27r + arcsin7. In the case of 
the microresonators (d > 1), the stable pinned fluxons have a "dip" before and near the inhomogeneity 
and the values in the dip are between arcsin7/(i and arcsin7, i.e., below the left asymptotic state. 

To complement and illustrate the analytical results in the previous sections, we have numerically 
solved the stationary equation (4) for the pinned fluxons and the corresponding linear eigenvalue 
problem (6) using a simple finite difference method and presented the results in Figures 23-26. Without 
loss of generality as far as stability is concerned, we depict the eigenvalues for o: = 0, i.e., A = A^. Thus 
an instability is indicated by the presence of a pair of eigenvalues with non-zero real parts. 

First, we consider the case of inhomogeneous Josephson junctions for a microresistor with d = 0. 
As is shown in Figure 17, when 7 = 0.15 and the defect length parameter L = 4.2, there are four 
possible pinned fluxons. In Figure 23, the numerically obtained profiles of pinned fluxons are a shown; 
all of them are clearly non- monotonous. The insets show the eigenvalues of the fluxons in the complex 
plane. Only one of them has no eigenvalues with non-zero real parts, confirming that there is exactly 
one stable pinned fluxon, which is non- monotonous for these parameter values. The four pinned fluxons 
belong to two different families, the ones with the smallest bump, i.e. panel (a), are on the blue curve 
and the others, i.e. panel (b), on the green curve in Figure 17. 




X X 

(a) (b) 

Figure 23: The four pinned fluxons admitted by the Josephson system with d = L = 4.2, and 
7 = 0.15. The insets show the eigenvalues of each fluxon; the top inset is related to the upper fluxon 
and the bottom inset to the lower fluxon. The vertical dashed lines show the edges of the defect. 

In Figure 17, the existence and the stability of the pinned fluxons for fixed d and 7 are presented 
in the (/i, L)-plane and it is shown that each pair of the fluxons collide in a saddle-node bifurcation at 
a critical L for a fixed 7. To complement these results, we take L = 4.2 and numerically follow the 
largest eigenvalue A = of the various fluxons when the induced current 7 changes. The results are 
shown in Figure 24. As before, the colouring corresponds to the colouring in Figure 17. Figure 24 
shows that there is a critical current for the existence of a pinned fluxon for a given length and depth of 
the inhomogeneity. The blue and green fluxons disappear in a saddle-node bifurcation. This happens 
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at a smaller value of 7 for the green fluxons (solutions in panel (b) in Figure 23) than for the blue 
fluxons (panel (a)). A physical interpretation of the saddle- node bifurcation is that the inhomogeneity 
is too short or long to pin a fluxon when the applied current exceeds the critical value. For 7 = 0.15, 
there are no red fluxons at this length, but they will exist for smaller values of 7. The red fluxons 
disappear when the fluxon "splits" in a homoclinic connection to 27r + arcsin7 and a blue pinned fluxon. 
Only one curve of red fluxons is visible. In theory, there is a second curve, but this exist in a tiny 
7-interval only and hence is not visible. 



Figure 24: The largest eigenvalue A = of the various fluxons as function of the induced current 7. 
The maximal eigenvalue at 7 = 0.15 of the fluxons in Figure 23 is at the intersection between the 
curves and the vertical dashed line. Note that the blue and green fluxons disappear in a saddle-node 
bifurcation. 

In Figures 25-26, we consider the case of a microresonator with d = 4. From the middle panel in 
Figure 22, it follows that there exist five pinned fluxons when 7 = 0.2, and L — 0.75. In Figure 25 
we show the numerically computed profiles of those pinned fluxons and their eigenvalues, where the 
colouring is as in Figure 22. The blue non-monotonic fluxon is stable while the blue monotonic one 
and red one are unstable. This confirms our analytical findings (see Theorem 11: there is at least one 
stable pinned fluxon). Moreover, it shows that there can be more than one stable fluxon: one of the 
fluxons on the black curve is stable too. So for d > 1, there is bi-stability for some values of L and 7. 

In Figure 26 we also present the critical eigenvalues of the five fluxons as a function of 7 when 
L = 0.75 is fixed. Similarly as in Figure 24, the pairs of blue and black fluxons collide in a saddle-node 
bifurcation, while the red fluxon breaks up at the maximal value of 7. 

For future research, it is of interest to expand our study to the case of two-dimensional Josephson 
junction with inhomogeneities. A particular example is the so-called window Josephson junction, which 
is a rectangular junction surrounded by an inhomogeneous 'idle' region with d = 0. The interested 
reader is referred to [4, 5, 8] and references therein for reviews of theoretical and experimental results 
on window Josephson junctions. Recently, fluxon scatterings in a 2D setup in the presence of a non-zero 
defect has been considered as well in [24]. 

One can also apply our method to study the existence of trapped solitons by inhomogeneities in 
Schrodinger equations, such as pinned optical solitons in a nonlinear Bragg media with a finite-size 
inhomogeneity (see, e.g., [13] and references therein) and trapped Bose-Einstein condensates by a finite 
square-well potential (see, e.g., [9, 23]). In general, the ideas presented in this paper are applicable 
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(a) (b) (c) 

Figure 25: The same as in Figure 23, but for a microresonator with (i = 4, L = 0.75, and 7 = 0.2, 
where there are five pinned fluxons. 




7 

Figure 26: The same as in Figure 24, but for the five pinned fluxons in Figure 25. 

to any system with locally (piecewise constant) varying parameters in the equations as can be seen in 
preprints by some of us [21] and [14]. 

Finally, the simulations in section 2 show how inhomogeneities can capture travelling fluxons. This 
suggests that the pinned fluxons analysed in this paper can be attractive or repelling, just as observed 
in [22] in case of the localised inhomogeneities. We are currently investigating the attractive and 
repelling interaction of the travelling fluxons with the pinned fluxons and will report on this in a future 
paper. 

Acknowledgement We would like to thank Daniele Avitabile for the use of his suite of simulation 
codes. 

A Largest eigenvalue of linearisation with no induced current 

Proof of Lemma 4 Let 7 = 0, d = 0, and fix the length L of the inhomogeneity. We denote the 
unique pinned fluxon with length L by (/)pin(^) (suppressing all other parameters). From (12), we see 
that (/)pin equals the sine-Gordon fluxon outside the inhomogeneity {\x\ > L) and the linearisation about 
the sine-Gordon fluxon is well-studied. The shifted pinned fluxon (j)pin{x) — tt is an odd function, hence 
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a quick inspection shows that the operator Cpin(x) is even in x (we suppress all other parameters in 
>Cpin)- All eigenvalues of Cpi^ are simple, thus Cpin{x) being even implies that all eigenfunctions are odd 
or even. The eigenfunction for the largest eigenvalue does not have any zeroes, thus this eigenfunction 
is even. 

For fixed A, the hnear ODE associated with (>Cpin — A)^^ = has two hnearly independent solutions. 
The asymptotic limits of (j)pin for x ±oc correspond to saddle points in the ODE (4) and the decay 
rate to these fixed points is like e~^. This implies that for A > —1, there is one solution of the ODE 
(>Cpin — A)^^ = that is exponentially decaying at +oc and there is one solution that is exponentially 
decaying at — oc. We denote the exponentially decaying function at — oo by V-(x; L, A). 

In [20], the linearisation of the sine-Gordon equation about the the fluxon 00 is studied in great 
detail. Using the results in this paper, we can derive an explicit expression for the solutions V-(x; L, A) 
(see also [10]), they are 

V-{x;L,0) = sech(x + x*(L)), x < -L 

i;_(x;L,A) = e^(^+^*(^)) [tanh(x + x*(L)) - /x], x < -L where = ^ATT, 

where x*(L) is given in Lemma 1. In the inhomogeneity, the linearised operator is simply Cpi^i — Dxx, 
hence the even solutions of Cpin — A are 

'^inhom(^; A) = Acos(\/— Ax), |x| < L, if A < 0; 

'^inhom(^5 0) — ^i \^\ < 

'^inhom(^; A) = A cosh(\/Aa:^) , \x\ < L, if A > 0. 

To have a continuously differentiable solution of (>Cpin — A)'0 = in i7^(R), we have to match V- and 
'^inhom ^ud its derivatives at x = — L (the conditions for x = L following immediately from this as the 
eigenfunction is even). This gives: 

• If A = (thus 1): 

A = sech^* and = — sech^* tanh ^* 

with ^* = —L + x*(L). This implies that ^* = and A = 0. From the relation for x*(L) in 
Lemma 1, it follows ^* 7^ only if L = 0, hence when there is no inhomogeneity. This confirms 
that the stationary sine-Gordon fluxon (the pinned fluxon for L = 0) has an eigenvalue zero, but 
none of the pinned fluxons with L > will have an eigenvalue zero for its linearisation >Cpin. 

• If A > (thus /i > 1), with = Ly^/x2 i ^nd again ^* = -L + x*(L): 

A cosh = e^^* [tanh^* - jj] 
-v/l -/xMsinhy* = e^^* [//(tanh^* - /x) + sech^^*] 

Hence /i (thus A) is determined by 



/i [tanh^* — /i] + sech^^* = — — 1 [tanh^* — jj] tanh^*. 

Using Lemma 1, this can be written as a relation between fi and (/)in (and hence fi and L as there 
is a bijection between (/)in G (0, tt) and L > 0): 



-/i 



fJ.-\-l ^2(1 + COS(/)in) + ^ (1 - cos (/)in) = 

r 1 / \ (22) 

V/i'-l W2(l + cos0in) tanh ""^^^ J. 
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It can be seen immediately that the right-hand side of (22) is positive. The left-hand side of (22) 
is always negative for > 1 as 

-/i (/X + T) + 1 - = 1 - /x^ - /iT - < -T - < 0, 



where we wrote T = ^ ^2(1 + cos(/)in), hence ^ (1 — cos 0in) = 1 — T^. Thus (22) has no solutions 
and there do no exist any eigenvalues A > 0. 



If -1 < A < (thus < fi <1), again with ^* = -L + x*(L) and now ?/* = L^/l - fi^: 



A cosy* = e^^* [tanh^* - ii] and \/l - /x^^giny* = e^^* [//(tanh^* - ^) + sech^^*]. 
Hence /x (thus also A) is determined by 

/i [tanh^* — /i] + sech^^* = \/l — ji^ [tanh^* — ji] tany*. 
Using the same relations as before, this can be written as a relation between ji and 



II +\ ^2(1 + COS(/)in) + ^ (1 - cos (/)in) = 

U + ^ V2(l + cos(/)in)l tan (^l^^^^^ 
Bringing all terms to the left and writing T(L) = ^ ^2(1 + cos 0in(i^)) G (0, 1) gives on the left 



F{L,ii) := -/X [/i + T] + 1 - + ^l-^i^ [^i + T] tan ( ^T^L 



Taking = 1 in this expression gives F(L, 1) = -T - < 0. If L < f , then F(L,0) = 1 - 

T^ + TtanL > as T G (0,1). If L > f , then c^-,^ < f and T > thus F(L, V^'-^ -gzi)!) > 

—2 + ^^^^2L^ tan(| — 5) = ^^(xe)' -^^^ ^ ^ 0. As L is fixed, we can choose e such that this 
expression is positive. Thus we can conclude that for all L > 0, there is at least one G (0, 1) 
that solves F(L, n) = 0. If L gets very large, then there will be many solutions, but we are 
interested in the largest one. 

□ 

Lemma 12 For 7 = and d — \ ^ e with £ small, the linearisation >Cpin(x; L, 0, 1 + s) about the 
monotone pinned fluxon (l)pin{x; L, 0, 1 + s) has a largest eigenvalue of the form sAi + 0{e^) with 



Ai 



sech^L [-L2sech^L(l + tanh^ L) + 2L tanhL(sech^L + 2(1 + sech^L)) + tanh^ L(6 + sech^L)] 



16(Lsech^L + tanhL) 

See Figure 27 for a sketch of Ai. Furthermore, if there are any other eigenvalues, then they must 
be near —1. Thus for 7 = and d close to 1, the monotone pinned fiuxons with d > 1 are linearly 
unstable. The nonlinear stability of Theorem 9 is confirmed by the sign of Ai for d < 1. 



Proof The monotone pinned fluxon for 7 = and d = 1 + £ with 6 <C 1 can be written as 

(/)o(x + £xl{L, 5)), X < -L, 



6pin(x;L,0, 1 + 5) 



0o(x) + (x; L, 5), |x| < L, 
(f)o{x — exl{L, 5)), X > L. 
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Figure 27: The eigenvalue factor Ai as function of L. 



Here (j)i{x]L^e) is an odd function satisfying 
^Dxx(j)i - (1 + 6) sin((/)o + e(j)i) + sin^o 
and x\[L^e) is such that (j)o(—L + exl) = (/)in 



0, Ixl < L, 



(23) 



(/)(){— L) + (—L). To find an approximation for 



and x^, we introduce the notation (f)^ = 00 (—L), and 0^ = 01 (—L), thus e(f)l 
length L is 



The half 



... a/2(1-: 



y^2(/i + (1 + e){l - COS(/>) Y^2(l - COS(/)in + (1 + £)(cOS(/)in - COS(/>)) ' 

where we used that h = e(cos(/>in — !)• With (/>in = (/>o + this becomes 



L 



L 



cos ( 



COS ( 



_e r 



(2(l-COS</>))3/2 



COS( 



(2(l-COS</>))3/2 



COS( 



+ 0{e^). 



COS( 



Rearranging this expression and using that cos 0o 
sechL 



1 — 2sech^(L), we get an approximation for 
[2L(l + tanh2L) -2tanhL] +0{e). 

Furthermore, x\ is given by (j)o{—L + ex\) = 0q + e(j)\. An expansion of (j)o{—L + ex\) shows that 

ssechL 



bl^exld^'^{-L) 



[2L(1 + tanh2L) - 2tanhL] -^0{e^). 



With 00 (—L) = 2sech(L), this shows that 



X 



^ = [2L(l + tanh2L) -2tanhL] +0(5). 



Next, we derive an approximation for the function 0i, using the differential equation (23). Expand- 
ing (23) in e gives 



D 



i)icos0o — sin 00 = 0{e) or >Co0i = sin0o + 0{e), 



(24) 
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with £0 = Dxx — cos 00. The homogeneous problem >Co'0 = has two independent solutions: '^^(x) = 
sechx and tpu{x) = xsechx + sinhx. In this, tpb{x) = |^0o(^) is bounded and ipu{x) unbounded as 
X ±00. By the variation-of-constants method, we find the general solution to (24), 

01 (x) = X sechx + A sechx B [x sechx + sinh x] 0(e), 

with A, S E M. As 01 must be odd, it follows that A = 0. Furthermore, the boundary condition at 
X = —L gives 0^ = —B (LsechL + sinhL) — LsechL + O(^), hence 

^ sechL(L tanh^ L — tanhL — 3L) 
4(L sechL + sinh L) 

Altogether we can conclude that 0i(x) = 0ii(x) + 0{s) with 

/ / \ 1 . sechL(L tanh^ L— tanh L— 3L) r i i • i i 
011 (x) = X Sechx + 4(LsechL+sinhL) ' ^^""^^ + Sl^"^ ' 

To find the largest eigenvalue of >Cpin(x; /i, 0, 1 + we will use perturbation theory. First we observe 
that for any L > 0, the linearisation Cq := >Cpin(x;L,0, 1) about the fluxon 0o has largest eigenvalue 
A = with eigenfunction is 0o. We have for \x\ < L 

Cpin{x; /i, 0, 1 + = Dxx - (l + e) cos(0o + £:0i) = jCo{x) - e (cos 0o - 0i sin 0o) + 0{e'^) 

and for x < —L 

>Cpin(x; /i, 0, 1 + 5) = >Co(^ + sxl) = >Co(x) + £:xi0o(x) sin0o + 0(^^). 
Thus the largest eigenvalue for >Cpin(x; /i, 0, 1 + 5) is A = + 6:Ai + 0{e^) and the eigenfunction is 



t/^ = 0Q + £'01 + 0{e ). The equation for Ai and '0i is 



£0^1 = Ai0o + /o(x), where /o(x) = < 



-X* sin0o(0'o)^ 



(cos 00 - 011 sin 0o) 00, 
x^ sin0o(^o)^ 
From (24) and the fact that >Co0o — O5 follows that 

>Co0ii = sin 00, hence >Co0ii = 2 (cos 0o — 0ii sin 0o) 0o 
>Co0o — O5 hence >Co0o — — sin( 

Thus, 



X < -L 
\x\ < L 
X > L 



(25) 



/o(x) Co < 



x^0o(x), X < -L, 
^0'ii(x), |x| < L, 
— x^0o(x), X > L. 



To find Ai, we multiply the eigenvalue equation (25) with 0q, integrate it, use integration by parts and 
>Co0o = and get 

POO 

Ai / {<P',f dx = 2x1 [i<PoiL)f - <i>'^iL)<i>'o{L)] - <P'li{L)<P',{L) + <^'n(L)</.'o'(L). 



with the explicit expressions for 0o and 0i, we get the expression in the Lemma. 

As the linearisation Cq about the sine-Gordon fluxon has exactly one eigenvalue (the one at zero), 
it follows immediately that if the perturbed linear operator has more eigenvalues, they have come out 
of the continuous spectrum, hence they are near —1. □ 
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